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ABSTRACT 

We present models of giant planet migration in evolving protoplanetary discs. We show that 
disc clearing by EUV photoevaporation can have a strong effect on the distribution of giant 
planet semi-major axes. During disc clearing planet migration is slowed or accelerated in the 
region where photoevaporation opens a gap in the disc, resulting in "deserts" where few gi- 
ant planets are found and corresponding "pile-ups" at smaller and larger radii. However, the 
precise locations and sizes of these features are strong functions of the efficiency of planetary 
accretion, and therefore also strongly dependent on planet mass. We suggest that photoevapo- 
rative disc clearing may be responsible for the pile-up of ~ Jupiter-mass planets at ~ 1 AU seen 
in exoplanet surveys, and show that observations of the distribution of exoplanet semi-major 
axes can be used to test models of both planet migration and disc clearing. 

Key words: accretion, accretion discs - planetary systems - planets and satellites: formation 
- protoplanetary discs - planet-disc interactions 



1 INTRODUCTION 

The study of planets dates back hundreds of years, but research 
in this area has accelerated dramatically in the decade and a half 
since the discovery of the fir st extra-solar planets (Ma yor & Ouelozl 
ll995tlMarcv & ButlerJ 199(3) . Ground-based survey s have now dis- 
cover ed over 500 planets around other stars (e.g.. lUdrv & Santosl 
2007), and the recent results from the Kepler satellit e have pushed 
our c ensus of exoplanets to over 2000 objects fe.g.. lBorucki et aU 
l201ll) . Many of the early exoplanets were found to orbit very close 
to their host stars, and it was immediately recognised that these "hot 
Jupiters" could not have formed at their current locations. Instead it 
is believed that planets form far from their hosts, at > 5-10AU, and 
migrate in wards due to tidal interactions with the protoplanetar y 
disc (e.g., iGoldreich & Tremaind [l980l : iLin & Papaloizoulll986l) . 
The migration of low-mass planets remains controversial, but so- 
called Type II migration, which applies to giant planets with masses 
> 0.5M Jllp , is now relatively well understood. Numerical models of 
Type II migration s how good agreement with the predictions of an- 
alytic theory (e.g.. iTakeuchi et alj|l99q) . and although the origin 
of hot Jupiters with periods of a few d ays remains uncertain (e.g., 
iRice et alj|2003 : IWu & Lithwic3l201lh . the observed distribution 
of giant planets with ~AU semi-major axes is b roadly consistent 
with the predictions of mig ration models (e.g., I Armitagel 120071 ; 
lAlexander & Armitagei2"0091) . 

How planet migration is stopped, however, remains an open 
question. The migration of low-mass planets in turbulent protoplan- 



* email: richard.alexander@leicester.ac.uk 



etary discs is subject to stochastic variations ^Nelson & Papaloizoul 
2004), and can be halted or even re versed in certain circumstances 
( Paardekooper & Papaloizou 2009). Type II migration, by contrast, 
continues as long as the disc accretes, and away from the inner 
disc edge it can only plausibly be halted by dispersal of the disc 
gas. Disc clearing at late times is thought to be dominated by 
photoevaporative winds, driven by ultraviolet and/or X-ray heating 
by the central sta r (e.g.. lAlexander et allboOfj ; iGorti et alj|2009l ; 
I Owen etal. 2010). These different winds disperse the gas disc and 
influence planet migration in quantitatively distinct ways, but ob- 
servationally it is not yet clear which he ating mechanism dominates 
jPascucci et al.l I20TTI ; ISacco et alj|2012h . Moreover, as our census 
of exoplanets has grown we have discovered strong features in what 
were previously thou ght to be sm ooth distributions of exoplanet 
properties (e.g.. IWright et al.i l2009j). In this Letter we use numeri- 
cal models of planet migration in evolving protoplanetary discs to 
investigate whether disc clearing leaves a characteristic signature 
on the observable properties of exoplanets. 



2 MODEL 

Our numerical model i s ess entially the same as that used in 
lAlexander & Armitage] d2009l . hereafter AA09); here we sum- 
marise only the salient details. We use a Monte Carlo approach, 
running large numbers of planet-disc models with randomly- 
sampled initial conditions to generate statistical distributions of 
planet properties. In each individual model the protoplanetary disc 
evolves due to viscous transport of angular momentum and photoe- 
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vaporation by the central star, while a single (giant) planet interacts 
tidally with the disc and undergoes Type II migration. The equa- 
tions governing the coupled evolution of the disc-planet system are 
(e.g jLin & Papaloizoulll986h 
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Here l.(R, t) is the disc surface density, t is time, R is (cylindrical) 
radius, v is the kinematic viscosity, and M, is the stellar mass. The 
first term on the right-hand side describes the viscous evolution 
of the disc, and t, w (R,t) is the mass-loss due to photoevaporation. 
The second term describes the response of the disc to the planetary 
torque A(R,a), which for a planet of mass M p = qM, a t radius 
(semi-major axis) a is given by (e.g jArmitage et all2002h 
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We adopt an alpha-prescription for the disc viscosity 

v(R) = aQH 2 , (5) 

where a = 0.01 is the Shakura & S unvaevl d 1973h viscosity pa- 
rameter and Q(R) = -\/GM,/R 3 . We choose a power-law form for 
the disc scale-height H oc R 5/ *, which gives a linear viscosity law 
v oc R. We normalize the power-law by setting the disc aspect ratio 
H/R = 0.0333 at R = 1 AU. 

As in AA09, we assume that photoevaporation is driven by 
extreme-ultraviolet (EUV) radiation from the central star, and 
make use of standard prescriptions for t,„(R). In our standard 
model we adopt a n ionizing flux O = 10 42 photons s (e.g., 
IPascucci & Ster zik 2009), which results in an integrated mass-loss 
rate of M w = 4x 10~ 10 M o yr~' for M, = 1M Q . When the inner disc 
is optically thick to ionizing photons (i.e., until a gap has opened in 
the disc), the mass-loss is concentrated around the critical radius 
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where R s = GM,/c 2 is the "gravitational r adius" and c s = 10km s 1 
is the sound sp eed of the ionized gas fHollenbach et al J 1 19941 ; 
iFont et alj2004l) . 

The initial disc surface density profile is given by 



2nR s R 



exp(-R/R s ), 



(7) 



where M d is the initial disc mass and the scale radius R s = 10AU 
effectively sets the viscous time-scale; for a = 0.01, this gives t v = 
R 2 /3v(R s ) a 5 x 10 4 yr. For each model Ma is drawn randomly from 
a log-normal distribution, with mean log 10 ((Md)/Mo) = -1.5 and 
a 3-cr spread of 0.5dex. In the absence of a planet the median disc 
lifetime is =4.5Myr, and AA09 showed that this underlying disc 
model is consistent with a wide range of observed properties of 
protoplanetary discs . Operationally, we solve Equation [TJ using a 
standard first-order explicit scheme, on a ^ 1/2 -spaced grid covering 
the range [0.04AU, 10000AU] (see AA09 for further details). 



We model accretion across the planetary orbit in terms of 
the efficiency parameter e, which is defined as the accretion rate 
on to the planet as a fr action of the disc accretion rate M disc (a) . 
Numerical studies (e.g.jLubow et alJI 19991 : Id'Angelo et alj|2002l; 
iLubow & d'Angelol [2006) have shown that e is a strong func- 
tion of planet mass, and we adopt the fitting formula derived by 
I Veras & Armitagd J2004h 
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The accretion rate on to the planet M p = e(M p )M l a x , and the 
accretion rate across the gap on to the inner disc is given by 
Mnner = M p /(l + e). Our results are not especially sensitive to 
the precise value of e max (we choose e max = 0.5), but do depend 
strongly on the form of e(M p ). We therefore also ran models with 
constant values of e(M p ), discussed below. 

We "form" giant planets by injecting a single planet of mass 
M p into each disc model at radius a p and time f p . We assume a 
constant formation radius a p , and the time of formation is assigned 
randomly within the range 0.25Myr < t p < t c (where t c is the time 
at which photoevaporation begins to clear the disc). We make no at- 
tempt to model the interaction between the migrating planet and the 
inner disc edge; instead, we simply remove planets which migrate 
to a < 0.15AU. As we wish to study how the properties of the sur- 
viving planets depend on their mass, we assign initial planet masses 
randomly in the range 0.5M ]l]p ^ M p ^ 5.0Mj up . This choice of a 
flat planet mass function serves merely to ensure we model a statis- 
tically significant number of massive planets, and is not intended to 
be representative of the observed exoplanet mass function (which 
declin es steeply with increasing planet mass; e.g., lUdrv & Santosl 
|2007|) . For each choice of input parameters we ran 1000 random re- 
alisations of the model. Typically around half of the planets migrate 
on to the star while the remainder survive, so each set of models re- 
sults in a statistical distribution of 400-700 planets. 

Our standard model considers a M, = 1M Q star and assumes 
that planets form at a p = 5AU. We also ran otherwise identical sets 



of models with constant e(M p )/e m 



1.0, 0.3, and a set of models 



with formation radius a p = 10AU. Finally, we considered a set of 
models with M, = 0.5M o . We again form planets (arbitrarily) at 
Op = 5AU, but in this case we re-scale the disc model by setting the 
mean disc mass log 10 ((M (i )/M o ) = -1.8, scale radius R s = 5AU, 
and ionizing flux O = 3 x 10 4 ' photon s~'. There is some observa- 
tiona l evidence for a linear scaling of dis c mass with stellar mass 
(e.g. jKlein et alj2003l:rScholz et alj2006h . but how the other model 
para meters scale with stellar mass in real systems is not well known 
(e.g jAlexander & Armitage1l2006l : lKennedv & Ken"vonll2009h . We 
have chosen these parameters, somewhat arbitrarily, to give a simi- 
lar median disc lifetime (== 4.5Myr) to the M, = 1M o case. 



3 RESULTS 

The radial distribution of surviving planets in the standard model 
is shown in Fig[TJ If we consider the complete sample of plan- 
ets we find no strong features in the distribution, and the overall 
shape of the di stribution reflects the average migration rate (e.g., 
lArmitagel2007h . For a given planet mass the Type II migration rate 
depends primarily on the disc viscosity and surface density, so the 
shape of the distribution in FigQ] primarily reflects our underlying 
disc model (and is consistent with AA09). The radial distribution of 
planets with final radii a < 4.5AU in the model with a p = 10AU is 
statistically indistinguishable from the standard model: this shows 
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Figure 1. Cumulative distribution of surviving planets in the standard 
model. We plot only the planets at a < 4.5AU, to eliminate artefacts in- 
troduced by the (arbitrary) choice of formation radius a p = 5AU. The thick 
black line shows the distribution for all planets, while the coloured lines 
denote sub-sets of this sample in different mass bins (as indicated in the 
legend). More massive planets show clear "deserts" in their radial distribu- 
tions, where few or no planets are found, close to the critical radius at which 
photoevaporation first opens a gap in the disc. 



that the distribution of planets in this region is determined primar- 
ily by migration, and is not very sensitive to our assumptions about 
where and when planets form. 

When we split the planets in the standard model into mass 
bins, however, several strong features appear. Most prominent are 
clear "deserts", where few or no planets of a given mass are found. 
The precise location of the desert is a strong function of planet 
mass: almost no planets are found between a ^ 1-2AU for the most 
massive planets (> 4M Jup ); between a =; 2-3. 5AU for 3-4M Jllp 
planets; and between a ^ 3-4AU for 2-3Mj up planets. There is 
weak evidence of a similar desert at even larger radii in the lowest 
mass bin (< 2M Jup ), but this is not statistically significant (and does 
not persist when the formation radius a p is increased). We also see 
significant excesses of planets ("pile-ups") both inside and outside 
these deserts, for all planet masses. 

These deserts and pile-ups are the result of the interaction be- 
tween the planet and the clearing gas disc, and can be readily un- 
derstood. As photoevaporation overcomes viscous accretion a gap 
opens in the disc at approximately R clit (i.e., at =i 1-2AU). Planets 
interior to this gap continue migrating for a short time while the 
inner disc drains, and consequently few planets are found just in- 
side ^ crit . Planets at larger radii suppress disc accretion, and if this 
suppression is strong enough the inner disc (inside R a i t ) becomes 
optically thin, allowing photoevaporation by direct irradiation to 
clear the outer disc. Planets which "trigger" disc clearing in this 
manner are subsequently unable to continue migrating, so we also 
see a deficit of planets just outside R CIlt . When integrated over many 
disc-planet models this results in a desert in the radial distribution 
of planets close to R alU and corresponding "pile-ups" at smaller 
and larger radii. The final location of any given planet, however, is 
a strong function of the rate of accretion across the planetary orbit 
(Mima). The photoevaporative wind is only able to open a gap in the 
disc once its rate exceeds the local disc accretion rate. If there is ef- 
ficient accretion of gas across its orbit (i.e., M mner ~ M disc , high val- 
ues of e) the planet has little effect on the gap-opening process, and 
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Figure 2. As Fig[T] but for models with constant e(M p ): e/e max = 1 (top) 
and e/e max = 0.3 (bottom). When gas is allowed to accrete freely across the 
gap (top) there are no deserts in the planet distributions, and no significant 
variations as a function of planet mass. By contrast, when accretion across 
the gap is suppressed (bottom) we see a clear lack of planets at a ss 1-3AU, 
and variations in the overall distribution with planet mass (because more 
massive planets migrate more slowly). 



consequently disc clearing has only a modest effect on migration. 
By contrast, low planetary accretion efficiencies severely limit the 
disc accretion rate interior to the planet's orbit (i.e., Mj nnel <k Mi\ x , 
small values of e), which accelerates the gap-opening process and 
has a much stronger effect on the planet's migration. We therefore 
see significant pile-ups of planets of all masses (though not always 
at the same radius), but find that deserts close to R ait are more pro- 
nounced for more massive planets. 

The sets of models with constant values of e(M p ) allow us to 
investigate this behaviour in more detail (see Fig0. For e/e max = 1 
(usually consistent with M p m 0.5M Jup ; see Equation[8]( there is co- 
pious accretion across the planetary orbit for all planets, and conse- 
quently disc clearing by photoevaporation does not leave a strong 
imprint on the final radial distribution of planets. Moreover, there 
are no significant differences in the radial distributions of planets 
of different masses: the distribution in each individual mass bin is 
statistically indistinguishable from the complete distribution. For 
e/emax = 0.3 (usually consistent with M p = 2.5Mj up ), however, 



©2012 RAS, MNRAS 000, [0f5] 



4 Alexander & Pascucci 



Exoplanet Data 



4 DISCUSSION 




a / AU 

Figure 3. Histogram of extra-solar planet semi-major axes. The shaded his- 
tograms show observational data for 131 single-planet systems with reflex 
velocities K > 20ms~' and host star masses in the range 0.75M o < M» < 
1.5M . The dotted histograms show the corresponding predictions of our 
standard model, normalised to match the numbers of planets in the observed 
samples. Although the agreement is not perfect, the most prominent feature 
- a significant excess of ~Jupiter-mass planets with semi-major axes be- 
tween 1-2AU - is clearly seen in both the model and the data. 



we see that almost no planets with M p > 2M Jup are found between 
1.5-3AU, and very few planets with M p < 2Mj„ p are found between 
0.6-2AU. There are significant pile-ups of planets inside these gaps 
for all planet masses. More significant pile-ups are seen for larger 
planet mass es: this is because mo re massive planets migrate more 
slowly (e.g., Syer & Clarke 1995), and are therefore more likely to 
survive at small radii (rather than migrating all the way on to the 
central star). These models clearly demonstrate that the crucial pa- 
rameter in determining the location of gaps in the radial distribution 
of planets is the accretion efficiency e, rather than the planet mass; 
the mass-dependence seen in the standard model (Fig[T) is in fact 
determined primarily by the form of e(M p ). 

In the models with M, = 0.5M o we again find large deserts in 
the radial distribution of massive planets, though the effect is only 
apparent for planets > 3M Jllp . The inner edges of the deserts are are 
at smaller radii than in the M, = 1M Q case, but we also see that the 
deserts are much wider, extending from a 1.5^1AU for 2-3Mj up 
planets, and from a ^ 1-3. 5AU for planets > 4M Jup . The smaller 
desert radii are expected from the linear scaling of ^ clit with stellar 
mass (Equation[6}, and represent a possible observational diagnos- 
tic. However, the wider deserts are primarily due to the lower disc 
masses in our model: massive planets open wider gaps in their par- 
ent gas discs, resulting in correspondingly wider deserts in the fi- 
nal distribution of planets. Without detailed knowledge of how disc 
evolution scales with stellar mass it is therefore difficult to make 
quantitative predictions for how the giant planet population varies 
around stars of different mass. Moreover, it is unclear whether this 
linear scaling with stellar mass can persist to lower-mass, M-type 
stars. For stars of mass M. ~ 0.1M G the critical radius R alt ~ 0.1- 
0.2AU, and is thus comparable to the radius of the inner disc edge. 
In this regime the manner in which photoevaporation affects disc 
clearing and planet migration is far from obvious. 



Our major result is that protoplanetary disc clearing by photoe- 
vaporation can have a strong influence on the distribution of giant 
planet semi-major axes. The general signature of this process is a 
gap or "desert" in the distribution of planets at radii close to the 
gap-opening radius R^t — 1-2AU, with corresponding "pile-ups" 
of planets immediately inside and outside this region. However, the 
precise locations and sizes of these features are strongly dependent 
on the efficiency of planetary accretion e, and consequently also 
on the planet mass. There is also a weak dependence on the un- 
derlying disc model, which determines the overall rate of planet 
migration. Numerical simulations have established that in the giant 
planet regime more massive planets accrete less efficiently, and our 
parametrization of e(M p ) (Equation [8} is base d on state-of-the-art 
numerical simulations of planetary accretion I Lubow et al.l 1 19991 : 
Id'Angelo et"al]|2002l : lEubow & d'Angeldl2006l . see also Fig.l in 
AA09), but it is clear from Fig[2]that the distribution of semi-major 
axes is rather sensitive to the exact form of e(M p ). We therefore pre- 
dict that the observed distribution of exoplanet semi-major axes can 
be used both as a diagnostic of disc dispersal, and as a test of the 
physics of planetary accretion. In particular, variations in the distri- 
bution as a function of exoplanet mass may be a useful means of de- 
termining how efficiently real planets accrete gas from their parent 
protoplanetary discs. We also note that the deserts and pile-ups we 
predict are the result of the complicated interaction between photo- 
evaporation, migration and accretion during disc clearing, which is 
usually neglected in simpler anal ytic calculations. More traditional 
popul ation synthesis models (e.g. Jlda & Lin 2008; Mordasi ni et al.l 
1 20091) typically adopt parametrized treatments of both migration 
and disc evolution, which necessarily neglect these more subtle ef- 
fects. Our results suggest that a more sophisticated approach may 
be required if such models are to make accurate predictions of ex- 
oplanet properties. 

Fig[3] shows a comparison between the predictions of 
our standard mo del and current e xoplanet data (taken from 
exoplanets . orq; IWright et alj|201 ll) . This compilation is not en- 
tirely free of biases, as it includes detections made using various 
different methods and selection criteria, but in the region of param- 
eter space which is most relevant here (a > O.IAU; M p > 0.5M Jup ) 
the overwhelming majority of the data come from radial velocity 
surveys. These data are essentially compl ete for reflex velocitie s 
K > 20 ms and periods P < 2000 days dCumming etai]|2008h . 
which for solar-mass stars roughly corresponds to a < 3AU and 
M p sini > 1.25M Jllp . The observed exoplanet sample is therefore 
not yet ideal for rigorous statistical testing of our models (as it is 
incomplete for masses < lMj up at ~AU radii), but is appropriate 
for qualitative comparisons between models and data. The major 
feature seen in the data in Fig[3]is a significant excess of -Jupiter- 
mass planets with semi-major axes between 1-2AU. This pile-up 
of planets is statist ically significant in single-planet systems (e.g., 
Wright et al. 200j^), and is not seen for higher planet massefl A 
similar feature is seen in our standard model: there is a significant 
pile-up of planets at 1-2AU, just interior to _R crit , which is only seen 
for planets < 3— 4Mj up (see also FigQ}. Although the agreement be- 
tween the data and our models is not perfect it is very encouraging, 



Strictly we are comparing the observed M p sin values with our true 
masses M p . However, the radii of pile-ups and deserts in our models depend 
primarily on e rather than M p , so the differences between M p and M p sini 
are essentially degenerate with the uncertainties in the form of e(M p ). 
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especially when one considers the large number of "un-tuned" free 
parameters in the model. 

An important simplification is our choice of photoevaporation 
model: we assume that the photoevaporative wind is driven by ion- 
izing extreme-ultraviolet (EUV) photons from the central star, with 
constant luminosity O. As R clit does not depend on O the assumed 
ionizing luminosity has no significant effect on planet migration, 
but some recent models suggest that photoevaporative disc clear- 
ing may in fact be primarily driven by X-rays and/or non-ionizing 
far-ul traviolet (FUV) radiation (e.g jGorti et al.f 2009; Ow en et al I 
l20ld) . These models also predict significantly higher photoevapo- 
ration rates, up to M w < 10~ 8 M o yr~'. Observationally it is not yet 
clear which is the primary heating mechanism, though recent stud- 
ies of diagnostic gas emi ssion lines do not favour models wher e 
X-ray heating dominates jPascucci et al.|[2oTTl ; ISacco et alj|2012h . 
Our reason for not including FUV or X-ray heating is simple. In 
both cases the photoevaporative flow is somewhat cooler than in 
the EUV case, leading to a larger critical radius R cnt . X-ray- and 
FUV-heated winds therefore tend to be launched from correspond- 
ingly larger radii, and their mass-loss profiles typically peak at 2- 
10AU or beyond. Unfortunately, this is close to the region where 
we be lieve giant planets form by core accretion (e.g.. |Pollack et al] 
1996), and consequently the effects of these winds on planet mi- 
gration are degenerate with our assumptions about where and when 
planets form. The general effect of X-ray and FUV photoevapora- 
tion is to create deserts in the distribution of planets between ~ 2- 
20AU, but without a fuller understanding of planet formation we 
cannot make detailed predictions about this process. We note in 
passing, however, that large photoevaporation rates (M w ~ 10~ 9 - 
10~ 8 M o yr~') imply much larger disc surface densities at the point 
when disc clearing begins than are found in our EUV-only mod- 
els. In this regime planet migration is rarely suppressed, leading to 
much faster planet migration and making it much less likely that 
planets will survive interior to R clit . Further work in this area is still 
needed, but our models suggest that very few giant planets should 
survive at small radii if photoevaporative clearing occurs at rates 
M w > lO^Moyr 1 . 

Finally, we note that photoevaporative clearing is not the only 
mechanism which can lead to features in the exoplanet distribu- 
tion. Various other effects, such as the location of the snow-line, 
formation in dead zones, or migration into planet "traps" (e.g ., 
iMasset et al]|2006l ; iMatsumura & Pudritdl2006l llda & Linl 2008), 
can also lead to pile-ups at specific radii. However, these features 
tend not to remain in fixed locations over the entire disc lifetime, 
and it is also not obvious if these mechanisms can reproduce the 
mass-dependence seen in the observations. We predict that proto- 
planetary disc clearing due to photoevaporation results in a pile- 
up of ~Jupiter-mass planets at ~ 1-2AU, and we suggest that 
this mechanism may be responsible for the similar pile-up seen 
in recent exoplanet surveys. As our sample of exoplanets becomes 
larger and more complete, their properties will provide important 
observational tests of the physics of disc clearing, planet migration 
and planetary accretion. 
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